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We investigate the electrostatic effects in doped topological insulators by developing a self con- 
sistent scheme for an interacting tight binding model. The presence of bulk carriers, in addition to 
surface electrons, generates an intrinsic inhomogeneous charge density in the vicinity of the surface 
and, as a result, band bending effects are present. We find that electron doping and hole doping 
produce band bending effects of similar magnitude and opposite signs. The presence of additional 
surface dopants breaks this approximate electron-hole symmetry and dramatically affects the mag- 
nitude of the band bending. Applying a gate potential can generate a depletion zone characterized 
by a vanishing carrier density. We find that the density profile in the transition zone between the 
depleted region and the bulk is independent of the applied potential. In thin films the electrostatic 
effects are strongly dependent on the carrier charge density. In addition, we find that substrate in- 
duced potentials can generate a Rashba type spin-orbit coupling in ultra thin topological insulator 
films. We calculate the profiles of bulk and surface states in topological insulator films and identify 
the conditions corresponding to both types of states being localized within the same region in space. 



I. INTRODUCTION 

The existence of metallic surface states with Dirac- 
hke dispersion represents the hallmark of topological 
insulators. In three dimensional (3D) topological insu- 
lators (TIs), the partially occupied surface states form a 
helical metal that has been predicted to host a wide range 
of new physical phenomena.'^^'^ Signatures of character- 
istic TI surface states have been observed in a family of 
strongly spin-orbit interacting Bi-hased materials.li2tSl 
However, most of the systems studied experimentally are 
not three dimensional topologically orderecP^ bulk solids, 
but rather doped TIs. Realizing bulk insulating samples, 
as well as separating the surface and bulk contributions 
in various types of measurements, represent serious chal- 
lenges. 

In doped 3D topological insulators, surface carriers 
that occupy the gapless surface states coexist with bulk 
carriers. As a result, the charge density in the vicinity of 
the TI surface is intrinsically non-homogeneous. Conse- 
quently, electrostatic effects, including band bending in 
the vicinity of the surface, contribute significantly to the 
low-energy physics in these systems. Additional electro- 
static effects may be generated by external electric fields 
that are present in various experimental setups involving 
TIs. Examples include applied gate potentials that cre- 
ate a bulk depletion zone and expose the surface states, 
gate potentials that control the bulk carrier density in 
thin films, and substrate-induced potentials that modify 
the band structure in epitaxially-grown films. 

In this paper we present a systematic theoretical study 
of electrostatic effects in doped TIs using an effective 
four band tight binding model^^ for the Bi-family topo- 
logical insulators. The Coulomb interaction is included 
at the Hartree level using a self-consistent scheme. We 
focus on systems with a slab geometry with (111) sur- 
faces. This surface orientation is the most relevant for 
experiment and is characterized by topological surface 
states with a localization length scale of the order of 



Inm. Note that other surface orientations are charac- 
terized by surface states with much larger localization 
length scales and, consequently, will exhibit significantly 
different electrostatic effects. The slab geometry allows 
us to study the effects of uniformly distributed bulk and 
surface dopants, e.g., those originating from the presence 
of Se vacancies. Our calculations ensure a realistic treat- 
ment of the density profiles originating from topological 
surface states, which are essential for understanding bulk 
screening effects and cannot be described within a sim- 
ple Thomas-Fermi approximation. The specific problems 
that we address have directly measurable consequences 
and we explicitly emphasize the links between our find- 
ings and experiment. 

Structure of the paper and main results 

In section II we describe the tight-binding model that 
we use for characterizing the low-energy properties of the 
doped TIs and the self-consistent method that we imple- 
ment numerically to account for the effects of Coulomb 
interaction. Our main results are presented in section III 
(for bulk TIs) and section IV (for TI thin films). Section 
V contains our conclusions. 

In Section III we investigate the electrostatic effects 
that occur in the vicinity of the surface of a bulk TI or 
near the interface between a TI and a trivial insulator. 
First, we address the question concerning the accuracy of 
Thomas-Fermi approximation in low-doped TIs (section 
III A) . We find that the Thomas- Fermi approximation is 
consistent with the numerical self-consistent calculations 
in regions with bulk charge density, but fails to describe 
the depleted region that may occur near the surface, as 
well as surface charge density contributions. Next, we 
study band bending in the vicinity of TI surfaces (section 
III B). The presence of this effect in TIs has direct exper- 
imental consequences, e.g., it generates differences^ be- 
tween surface-sensitive probes, such as the angle-resolved 
photoemission spectroscopy, and bulk-sensitive probes. 
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for example Shubnikov-de Haas measurements. In doped 
TIs, band bending is the result of the bulk electrons being 
pushed away by the surface charge of the helical metal. 
In addition, surface dopants, such as Se vacancies^^"^^ 
that migrate to the surface, can significantly modify the 
bending of the bands at the boundary. We investigate the 
difference between band bending in electron-doped and 
hole-doped TIs, as well as the role of surface dopants. 
To facilitate visualizing the band bending effects, we cal- 
culate the local density of states (LDOS) in the vicin- 
ity of the surface. Band bending, together with the 
charge density contributions of surface states, generate 
intrinsically non-homogeneous density profiles and effec- 
tive electrostatic potentials that vary significantly in the 
vicinity of the system boundaries. These potentials can 
be probed experimentally, for example using second har- 
monic generationl ^^ l ^^l We also study the charge density 
profiles in the presence on an external field. We find that 
applying a gate potential can create a depletion zone near 
the surface characterized by a vanishing bulk carrier den- 
sity, while deep inside the bulk the potential is screened 
and the charge density is constant. We also find that 
the density profile in the transition zone between the de- 
pleted region and the bulk is independent of the applied 
potential. Creating a depleted region could allow access 
to the surface states without interference from the bulk 
carriers. 

In section IV we investigate different types of electro- 
static effects that can occur in topological insulator thin 
films. In general, screening in TI films is significantly 
weaker that in bulk systems. Consequently, the band 
bending effects can be quite severe and, typically, the 
system does not have a bulk region characterized by con- 
stant charge density. First, we study the effects of an 
external potential, which may occur, for example, as a 
result of a charge transfer from the substrate^ (section 
IV A) . We find that the substrate-induced band bending 
is strongly dependent on the carrier density. In addition, 
in ultra-thin TI films a gap opens at the Dirac point as a 
result of the overlap between the surface states from the 
two surfaces and, in the presence of a substrate-induced 
potential difference, an effective Rashba-type spin orbit 
splitting is generated.'^ Finally, we calculate the profiles 
of the surface and bulk states in TI films (section IV 
B) and show that, in certain conditions, the lowest en- 
ergy states from the conduction band become localized 
near the TI surface. This opens the interesting possibility 
of coupling surface and bulk states, which typically are 
spatially separated, using, for example, optically-induced 
transitions. 



trons occupying the topological surface states is screened 
by the bulk charges. To study the resulting charge den- 
sity profile we consider a minimal model described by the 
Hamiltonian: 
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The first term, H^i, is an effective tight-binding model 
on a triangular lattice with 4-bands per quintuple layer. 
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where c|„ = are spinors and creates 

an electron in one of the 4 bands labeled by a = 1,2 
and a =t,i- The first and the second term of Hti 
describe the intra and inter band hoping respectively. 
The tight binding model parameters used in this calcula- 
tion are the intra layer hoping elements t\^^ = 1.42995 
and i,-^-* = 0.0318196, the inter-layer hoping elements 

-0.0413289 and band edges 



t''/^ = -2.95299 and 
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17.6495 and the spin orbit cou- 



plings A(i) ^ 0.291864 and A(2) = 0.115223 all measured 
in eV. The Coulomb repulsion. Hint, is described at the 
mean field level by the term. 
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where = J2a,a fLaf^ao- is the local electron density 
and riion is the ionic density. 

In a slab geometry, the in-plane wavevector ky is a 
good quantum number, due to the translational invari- 
ance along the layers. The local density, (fij), depends 
only on the layer index and can be expressed as the av- 
erage of contributions coming from the different momen- 
tum sectors. For a slab of N layers, each sector contains 
AN orbitals and is diagonalized independently. For a 
given density profile, the mean field local chemical po- 
tential can be determined by solving self consistently the 
discrete one dimensional Poisson equation, 



fJ-i = ^J■ax + Sfii i{ni)) 



(4) 



where 



N 



Sfii ((77,)) = 2[/|| {Sni) + C/_L E I' - ^'1 (^^') (5) 



II. MODEL AND METHOD 

Topological insulators from the Bi family have a lay- 
ered structure consisting in quintuple layers of Bi and Se 
(or Te) atoms oriented parallel to the (111) plane. In a 
slab with (111) surfaces, the potential created by elec- 



vanishes at uniform electric density, {6ni) — (rii) — rig, rif. 
the total electron density and 

flo,l^^^o~El^P^l- -Si, (6) 



is a constant external potential, where the magnitude of 
the external electric field E, the parabolic potential P 



3 



and a surface potential Si are given by 

E = Eext + Uj_{5nsurf,l — Sllsurf.NL) C^) 

P ^ U^{n,on~n,) (8) 

Si = 2U\\[5l^i5nsurf.l+ Si^NL^nsurf^NL]- (9) 

Here Eext represents the external electric field, dusurf,!, 
Snsurf,N the doping of the first and last surface respec- 
tively measured relative to the ion density, and U±, 
[/|| are respectively the intra and inter layer Coulomb 
repulsion. The values of the coupling constants are 
Ui— — and = 7-C/i where a is the lattice con- 
stant of each layer, c the inter layer distance, 7 « 0.1936 
a geometrical factor and e sa 90 the electric permitivity 
of the bulk. This corresponds to a value of the interac- 
tion parameter J7|| = 4ey which is the one that we will 
use throughout this paper. 



Carrier density 

In this paper we will present a self-consistent mean- 
field calculation of screening effects in topological insula- 
tors that describes the density profiles generated by the 
contributions of both surface and bulk states. In a slab 
with iVi layers the carrier density is a sum of two contri- 
butions, the bulk doping and the surface doping, 

r r , ^''^surf.l + Snsurf,N /, „^ 

dne = dn^on-{ ^—77 —■ (10) 

where Sriion and Susurf is the bulk and surface doping 
respectively. The bulk doping typically comes from Se 
vacancies. Their volume density, 5pion, varies between 
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10 
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which translates in a per unit-cell dop- 



ing of Suion = Vu.c.Spion ~ 10 - 10 , whcrc Vu.c. ~ 
1.42 X lO^^^cm"^ is the volume of the unit cell. The 
surface doping comes from Se vacancies on the surface 
created during the cleaving process or donors that come 
from the solution in which the sample is stored. Their 
surface density Sasurf typically varies between 10 
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cm 



and lO^'^cm ^. This corresponds to a per-surface layer 



fVu., 



10--^ - 10- 



which can 



doping of Susurf = j/g- 

be 10 — 1000 times bigger than the bulk doping at the 
surface layer. 



III. ELECTROSTATIC EFFECTS NEAR THE 
SURFACE OF A BULK TOPOLOGICAL 
INSULATOR 
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Figure 1: (color online) The energy spectrum for p = 6 x 
10"* as a function of |k|||, featuring the occupied (red) and 
unoccupied (grey) states of the conduction band, the valence 
band states (green) and the Dirac cones (purple) which are 
separated due to the bias. 



total charge of the system. At the surface, the localized 
charge from the topological states repels the bulk charge 
and, as a result, a partly depleted zone zone separates the 
surface and bulk regions. In the surface region and in the 
depleted zone the electric field in nonzero. The presence 
of electric fields near the surface results in the bending of 
the bands upwards or downwards, depending on the sign 
of the bulk carriers and of the surface dopants. Under- 
standing the detailed band bending mechanism is essen- 
tial for the interpretation of various surface and bulk sen- 
sitive measurements. In this paper we perform a system- 
atic study of these effects this by solving the electrostatic 
problem self-consistently. To emphasize the necessity of 
a self consistent scheme that accounts for the interaction, 
we compare the density profiles obtained in the presence 
and in the absence of interactions. As shown in Fig. [2] 
for a 60 layer slab with bulk density of 8 x 10^^/cell, the 
presence of interactions results in the development of a 
depletion zone near the surface. This non-homogeneous 
density profile generates a position-dependent electro- 
static potential that modifies the spectral properties of 
the system. 



Topological insulator samples exhibit very rich elec- 
trostatic behavior, especially near the surface. This is 
mainly due to the presence of bulk impurities, e.g., Se 
vacancieJiSl^ which results in an excess charge in the bulk. 
The surface itself may contain charged impurities, due to 
air exposure or due to the migration of Se vacancies from 
the bullP^l^ which provides additional contributions to the 



A. Accuracy of Thomas-Fermi approximation 

In metals and semiconductors screening effects are ac- 
curately described in the context of the Thomas-Fermi 
(TF) approximation in which the local charge density 
p (x) at position x is a function of the local chemical 



4 




_8 I \ \ \ \ \ 1 

5 10 15 20 25 30 



Figure 2: (color online) The local excess charge density, 
pi — pBuik in the vicinity of the surface with and without elec- 
trostatic effects for a system with bulk density 8 x 10~^/cell. 
In the absence of electrostatic effects, the density of bulk car- 
riers in a charge neutral system almost uniform. Small varia- 
tions are present in a thin surface region consisting in a few 
quintuple layers. In the presence of electrostatic effects, the 
charge density associated with surface states pushes the bulk 
carriers away from the surface creating a depletion zone. As a 
result, the electric field becomes nonzero in wide region that 
extends more than 20 quintuple layers from the surface. 

potential fi (x) 

p(x)« \,ondie)de, (11) 

where eb is the band edge and Ncond (e) the density of 
states of the conduction band. The electrostatic problem 
is then reduced in solving the Poisson equation for /i (x) . 
This approximation is violated at low densities, when the 
screening length ^ is much smaller than the de-Broglie 
wave length of the conduction electrons. For topological 
insulators this happens near the surface, where the topo- 
logical surface charges repel the bulk charges to create 
a charge depletion zone as shown in Fig. |2] Further- 
more, the profiles of the surface states extend up to 3-4 
layers into the bulk and need to be determined explic- 
itly by solving a quantum problem. To verify the va- 
lidity of the Thomas-Fermi approximation, we calculate 
the density profiles and the local electrochemical poten- 
tials for a thick TI slab in the presence of an external 
bias potential of varying strength and compare the re- 
lation between these two quantities with that predicted 
by the Thomas- Fermi approximation. The results are 
shown in the left panel of Fig. [3j At high densities, the 
self consistent p vs /i curve agrees well with the Thomas- 
Fermi approximation. However, at low densities (less 
than 5 x 10~^/cell) the self-consistent scheme predicts a 
small residual charge in the depletion zone that is absent 
in the Thomas-Fermi approximation. To shed light on 
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Figure 3: (color online) Left: The self consistent local charge 
density pi plotted against the local /i; , for a 60 layer slab with 
various densities and external biases /S.V = — OAeV . The 
solid line is the Thomas- Fermi approximation. The discrep- 
ancy is larger at low densities where there is a residual charge 
density due to valence band contributions. Right: Conduc- 
tion band charge density versus local potential. Note that 
the pconduction vs /I curves are in better agreement with the 
Thomas-Fermi approximation everywhere except the tail of 
the transition zone. 

the origin of the discrepancy, we show the self-consistent 
Pconduction VS fx relation for the conduction band elec- 
trons, without contributions from the valence band. The 
result, shown in the right panel of Fig. |3] indicates a bet- 
ter agreement with Thomas-Fermi at low densities. We 
conclude that, in addition to the surface contribution due 
to the characteristic helical metal, the Thomas-Fermi ap- 
proximation cannot capture the physics of the depletion 
zone, which is dominated by valence band contributions. 
In essence, this discrepancy is a quantum effect and is 
due to the change in the valence state profiles induced 
by the nonuniform effective potential near the surface. 
Taking this effect into account may be important for the 
understanding of surface phenomena and for interpreting 
surface-sensitive measurements. 



B. Band-bending 

Understanding band bending in doped TIs is critical 
for explaining the differences^^ between surface-sensitive 
probes, such as the angle-resolved photoemission spec- 
troscopy (ARPES), and bulk-sensitive probes, for exam- 
ple Shubnikov-de Haas measurements. In general, band 
bending in the vicinity of the surface is present in all 
doped TI s ystem s. However, surface dopants, such as 
Se vacancie^^^EH that migrate to the surface, can signif- 
icantly modify the bending of the bands at the bound- 
ary. Note that the concentration of carriers generated by 
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Figure 4: (color online) The local charge density for elec- 
tron doped and hole doped TIs in the vicinity of the surface. 
The system has no surface dopants. The corresponding band 
bending near the surface has opposite directions for electron 
doped and hole doped systems and there is an approximate 
symmetry between the two cases. The surface charge (corre- 
sponding to the first quintuple layer) is out of scale. 



surface Se vacancies can be controlled by depositing O2 
on the surface, which has been demonstrated^" to be an 
electron acceptor on Bi2Se3 surfaces. 

Electron vs hole-doped topological insulators. In the 
absence of surface dopants, the surface charge density 
profile exhibits an approximate electron-hole symmetry, 
as shown in Fig. |4] The depletion zone described 
above represents a region with lower carrier concentra- 
tion, which corresponds to a lower electron density in 
n-type doped TIs and a higher electron density in p-type 
systems (see Fig. |4]). These highly nonuniform density 
profiles can be understood in terms of band bending in 
the vicinity of the TI surface. Rigorously speaking, band 
bending can be associated with a position-dependent lo- 
cal density of states (LDOS). The LDOS in the vicinity 
of the TI surface for both electron doped and hole doped 
systems is shown in Fig. [5] The approximate electron- 
hole symmetry that characterizes the system in the ab- 
sence of surface dopants is manifest in in Fig. 



5(a) and 

|5(b)[ Quantitatively, in electron-doped systems with car- 
rier densities of the order of lO^^cm"'^ the band bending 
near the surface is positive and is characterized by an en- 
ergy shift of about QOmeV. The p-doped case is charac- 
terized by a negative band bending of similar magnitude. 

In the presence of surface dopants, the approximate 
electron-hole symmetry is broken. In particular, as 
shown in Fig. 5(c) and |5(d)[ positively charged surface 
dopants (i.e., donors) strongly push the direction of the 
band bending downwards. As seen in the left panel of 
the figure, an appropriate density of surface donors can 
basically eliminate band bending in an electron doped. 
By contrast, adding donors on the surface of a p-doped 
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Figure 5: (color online) Local density of states as a function 
of the distance from the surface in electron-doped (panels a 
and c) and hole-doped (panels b and d) TIs for 8 x 10~^ 
dopants/cell. The top panels (a and b) show the band bending 
in systems with no surface dopants, while the bottom panels 
(c and d) correspond to a system with 2.5 x 10"'' electron 
donor impurities on the surface. Note that the surface donors 
almost annihilate band bending in the n-doped system (panel 

c) and strongly enhances the effect in the p-doped TI (panel 

d) . Adding electron acceptors on the surface will have an 
opposite effect. 



TI will further enhance the downward band bending (see 
Fig. |5(d) I. Of course, adding electron acceptors on the 
surface will have opposite effects, i.e., will enhance band 
banding in electron-doped systems and reduce or even 
reverse the effect in hole-doped TIs. 

Density profiles in the vicinity of the surface of a topo- 
logical insulator. The intrinsically non-homogeneous 
density profiles induced by the surface charges gener- 
ate an effective electrostatic potential that varies sig- 
nificantly in the vicinity of the TI surface. A powerful 
method to probe this potential is the second harmonic 
generation (SHG) using ultrafast laser pulses.'^'^ Fur- 
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Figure 6: (color online) Charge density profiles for an 80 layer 
n-doped TI slab with 8 x 10^^ carriers per unit cell and no 
surface dopants. The different lines correspond to different 
gate voltage AV. Note that, increasing the system size will 
result in expanding the bulk region characterized by a con- 
stant density p = 0, but the density profiles in the surface and 
the depletion (excess) regions will not be affected. Applying a 
negative gate voltage results in the expansion of the depletion 
zone, while a strong positive potential can creates a region 
with excess carriers near the surface. The profile of the tran- 
sition zone (i.e., the region between the depletion region with 
no carriers and the bulk region with constant density) is in- 
dependent on the applied potential, as long as depletion/bulk 
regions exist, which is consistent with the Thomas-Fermi ap- 
proximation. 



thermore, applying gate voltages, which can allow the 
control of the charge carrier density in the vicinity of the 
surface, generates external fields that are nonzero in the 
depletion and surface regions. To determine the effect of 
a gate potential on the TI carrier density, we determine 
self-consistently the density profiles for different values of 
the gate potential. The results are shown in Fig. [6] The 
calculations were done on an 80 layer TI slab. However, 
if a bulk region characterized by constant charge density 
(i.e., a neutral region) exists, the density profiles outside 
the bulk region are independent on the size of the system. 
Hence, our result are relevant for bulk TI crystals with 
(111) surfaces. Applying a negative gate voltage results 
in the expansion of the depletion zone. However, we find 
that the profile of the transition zone that separates the 
depletion region with no carriers from the bulk does not 
depend on the applied potential. This is consistent with 
the predictions of Thomas-Fermi approximation, which is 
not surprising if we consider our findings concerning the 
validity of this approximation in the transition regime 
(see Fig. |3|. Finally, we note that a positive gate poten- 
tial reduces the depletion zone and, if strong enough, can 
even create a region with excess carriers near the surface. 

Next, we address the question concerning the depen- 



dence of the energy of the Dirac points and the depen- 
dence of the gap edges on the applied gate potential. The 
results for both n-doped and p-doped TIs, with and with- 
out surface dopants, are shown in Fig. |7] As expected, 
the effective gap is reduced by the external potential. The 
effect is also dependent on the concentration of surface 
dopants. Note that, depending on the type of bulk and 
surface dopants, as well as the applied gate potential, the 
Dirac points can be inside the bulk gap or buried inside 
either the valence or the conduction band. Using the 
right combination of dopants and applied gate potentials 
is key whenever exposing the Dirac point is critical for 
probing certain physical properties, e.g., the opening of 
a gap at the Dirac point. 



IV. ELECTROSTATIC EFFECTS IN 
TOPOLOGICAL INSULATOR THIN FILMS 

In thin TI films, i.e., films containing ^ 10 layers 
(lOnm), most of the carriers may come from surface 
donors. Consequently, it is possible that the "bulk" is 
no longer charge neutral. Strictly speaking, in most thin 
TI films there is no bulk region as defined above, i.e., a 
neutral region characterized by a uniform charge density. 
Also, external potentials can create strong perturbations, 
as screening is significantly weaker than in bulk TIs. All 
these characteristics lead to certain specific electrostatic 
effects that are not present in bulk systems. 



A. Substrate induced band bending in thin films 

Thin films are often grown on substrates using, for exam- 
ple, molecular beam epitaxy (MBE). In epitaxially grown 
films, a potential difference between the two surfaces of 
the TI may occur as a result of the charge transfer from 
the substrate.^ Furthermore, in the ultra-thin film limit 
(less than 6 quintuple layers) the overlap between the 
surface states from the two surfaces of the film cannot 
be neglected and, as a result, a gap opens at the Dirac 
point. In the presence of a substrate-induced potential 
difference, an effective Rashba type spin-orbit splitting is 
generated at k ^ 

First, we study a doped ultra thin TI film in the pres- 
ence of a substrate-induced potential difference. The 
spectra for films with six quintuple layers or less are 
shown in Fig. [8] Note that the Rashba-type band spit- 
ting and the gap opening in the surface spectrum due to 
the overlap of the surface states localized near the two 
surfaces are in good agreement with the experimental 
finding^lS. 

Next, we address the question of band bending in 
thin films. As mentioned above, in thin films the ex- 
ternal fields are weakly screened. This provides an ef- 
fective way to control the chemical potential in doped 
TI thin films and to realize bulk-insulating systems by 
applying gate voltages to remove the excess bulk charge 



7 



p-doped, no surf, dopants 

\l i ! 




n-doped, no surf, dopants 



0.2 - 



N =1 N =2 N =3 N =4 N =5 N =6 




> 



C 
LU 



0.2 


-0.2 
-0.4 




0.1 0.2 0.3 

AV(eV) 



0.1 0.2 0.3 

AV{eV) 



-0.2 0.2-0.2 0.2-0.2 0.2-0.2 0.2-0.2 0.2-0.2 0.2 
ka ka ka ka ka ka 

Figure 8: Ultra thin film spectra in the presence of a 
substrate-induced potential. Note the Rashba type band 
splitting splitting for A^^ > 3. Because the surface states 
extent about three layers away from the film surface, they 
hybridize with the surface states on the opposite side, which 
gives rise to a gap. 
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Figure 7: (color online) Energy of the Dirac points and gap 
edge energies as functions of the applied gate potential. To 
help the eye, we draw hues corresponding to: the bottom 
of the conduction band (blue), the top of the valence band 
(magenta) and the upper (red) and lower (green) Dirac pints. 
We note that, although the exact position of the Dirac cones 
relative to the bulk bands is model dependent, the general 
trends shown here are expected to hold in general. 



carriers .l^i^ We note that, in systems characterized by 
weak substrate-induced band bending, it is possible to 
have the Dirac points associated with the two surfaces 
approximately at the same energjEl. After removing the 
excess carriers, the film is no longer charge neutral and 
an effective potential with a parabolic layer dependence 
is created inside the film. Nonetheless, for typical film 
sizes and initial carrier concentrations [e.g. 10 layers. 



Figure 9: (color online) Local density of states for a 10 layer 
Tl film in the presence of an external potential difference, 
AV = O.leV. (a) Sn^ = 0, (b) Sn^ = lO-V(unitcell) (c) 
Srig = 5 X 10~^/(unitcell). Note that the band bending effect 
depends strongly on the carrier concentration. 



8 X 10^^ electrons/cell] the corresponding band bending 
effects are negligible. For comparison, we note that, for a 
TI film with an initial carrier concentration of 8 x 10^^, 
after completely depleting the excess charge the gap will 
close due to band bending effects if the thickness of the 
film exceeds 60 layers. 

In addition to the effective potential generated by the 
removal of charge carriers, external gate potentials and 
substrate-induced potential differences can significantly 
modify the spectrum. In Fig. |9] we show the self- 
consistently calculated band bending of a thin film with 
no surface dopants for three different values of the car- 
rier density. Note that, for a given potential difference, 
the intensity of the band bending effect depends strongly 
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cupied surface states. These surface and bulk-type states 
could be coupled optically, if the corresponding matrix el- 
ements are nonzero. Typically, the main problem in TIs 
is that the surface states are localized within a few quin- 
tuple layers from the surface while the bulk states are 
extended, hence the overlap is vanishingly small in the 
thermodynamic limit. However, in the presence of strong 
negative band bending, the low-energy conduction states 
become localized near the surface. This situation, which 
ensures the necessary condition for a strong surface-bulk 
overlap, is illustrated in Fig. [TO] Future studies are 
necessary to determine explicitly the optical matrix ele- 
ments associated with transitions from the surface states 
to the conduction band in doped TIs with strong band 
bending. 



Figure 10: Profiles of surface and bulk states for the same 

In this case the Dirac point is 



parameters as in Fig 5(d) 



buried inside the valence but the surface Fermi line is exposed 
due to band bending. The substantial spatial overlap may 
allow optical transitions between occupied surface states and 
empty conduction states. 



on the carrier density. As expected, higher carrier con- 
centrations provide better screening, which results in a 
weaker band bending. Also note that, depending on the 
carrier density, the chemical potential may be inside the 
gap (for films with a density close to half filling), it may 
cross the bottom of the conduction band at intermediate 
doping, or it may be completely inside the conduction 
band at large filling. 



B. Surface states and bulk states in topological 
insulator films 

In TIs the surface states associated with (111) surfaces 
are localized within about one quintuple layer from the 
surface and, consequently, the matrix elements of any 
operator between surface and bulk states vanish in the 
thermodynamic limit. However, in thin films the overlap 
between surface and bulk states becomes finite. More- 
over, we show that, in the presence of a local electric 
field that pushes the electrons toward the surface, the 
lowest energy bulk states become localized in the vicinity 
of the TI surface. This opens the possibility of optically 
coupling surface and bulk states in topological insulators 
with finely tuned surface and bulk density. One possibil- 
ity is to use a p-doped TI with strong downward band 
bending near the surface (See Fig |5(d)| . The tip of the 
Dirac cone will be buried inside the valence band, but 
the Fermi line corresponding to the surface states will 
be exposed due to the strong band bending. For any 
given wave vector in the vicinity of the surface Fermi k- 
vector, there are empty conduction states above the oc- 



V. SUMMARY AND CONCLUSIONS 

We use a self consistent scheme to study electrostatic 
effects in topological insulators. We find that in the vicin- 
ity of the TI surface a depletion zone characterized by 
the absence of bulk carriers may occur. We find that the 
Thomas-Fermi approximation fails in this depletion zone 
and in the surface region where the contribution from the 
topological states is dominant. As a result of the intrin- 
sic charge inhomogeneity in the vicinity of the surface 
we find that band bending effects are typically present. 
In the absence of surface dopants the band bending for 
electron and hole systems have similar magnitudes and 
opposite sign. The presence of surface dopants breaks 
this approximate electron-hole symmetry and dramati- 
cally affects the magnitude of the effect. We find that 
applying an external gate voltage can expand the deple- 
tion zone away from the surface. The transition region 
between the depletion zone and the bulk is character- 
ized by a charge profile that is independent of the ap- 
plied potential, in accordance with the Thomas-Fermi 
approximation. In thin films there is less bulk charge 
and, consequently weaker screening effects are observed. 
We show that band bending due to external fields de- 
pends strongly on the total charge density of the sample. 
We find that substrate induced potentials can generate 
a Rashba type spin-orbit coupling in ultra thin TI films. 
We calculate the surface and bulk state profiles and find 
that, in certain conditions, they are localized within the 
same region in space. Future work is needed determine 
if selection rules allow optical transitions between these 
states. 



Acknowledgments 

We would like to thank Dennis Drew for stimulating 
and useful discussions. 



9 



^ C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

^ B. A. Bernevig and S.-C. Zhang, Physical Review Letters 

96, 106802 (2006). 
^ L. Fu, C. L. Kane, and E. Mele, Phys. Rev. Lett. 98, 

106803 (2007). 

" L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
^ J. E. Moore and L. Balents, Physical Review B 75, 121306 
(2007). 

® R. Roy, Phys. Rev. B 79, 195321 (2009). 
X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 78, 
195424 (2008). 

® A. M. Essin, J. E. Moore, and D. Vanderbih, Phys. Rev. 

Lett. 102, 146805 (2009). 
^ L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 
^° B. Seradjeh, J. E. Moore, and M. Franz, Phys. Rev. Lett. 

103, 066402 (2009). 
" D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, 

and M. Z. Hasan, Nature 452, 970 (2008). 
^2 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Ban- 

sil, D. Grauer, Y. S. Hor, R. J. Cava, et al., Nat. Phys. 5, 

398 (2009). 

D. Hsieh, Y. Xia, D. Qian, L. Wray, F. Meier, J. H. Dil, 
J. Osterwalder, L. Patthcy, A. V. Fcdorov, H. Lin, et al., 
Phys. Rev. Lett. 103, 146401 (2009). 



X.-G. Wen, Int. J. Mod. Phys. B4, 239 (1990). 

J. A. Hutasoit and T. D. Stanescu, Phys. Rev. B 84, 

085103 (2011). 

1® J. G. Analytis, J.-H. Chu, Y. Chen, F. Corredor, R. D. 
McDonald, Z. X. Shen, and I. R. Fisher, Phys. Rev. B 81, 

205407 (2010). 

D. Hsieh, J. W. Mclver, D. H. Torchinsky, D. R. Gardner, 
Y. S. Lee, and N. Gedik, Phys. Rev. Lett. 106, 057401 
(2011). 

1® D. Hsieh, F. Mahmood, J. W. Mclver, D. R. Gardner, Y. S. 

Lee, and N. Gedik, Phys. Rev. Lett. 107, 077401 (2011). 

Y. Zhang, K. H. andCui Zu Chang, C.-L. Song, L.-L. 

Wang, X. Chen, J.-F. Jia, Z. Fang, X. Dai, W.-Y. Shan, 

S.-Q. Shen, et al., Nat. Phys. 6, 584 (2010). 
2° Y. L. Chen, J.-H. Chu, J. G. Analytis, Z. K. Liu, 

K. Igarashi, H.-H. Kuo, X. L. Qi, S. K. Mo, R. G. Moore, 

D. H. Lu, et al.. Science 329, 659 (2010). 

21 J. Chen, H. J. Qin, F. Yang, J. Liu, T. Guan, F. M. Qu, 
G. H. Zhang, J. R. Shi, X. C. Xie, C. L. Yang, et al., Phys. 
Rev. Lett. 105, 176602 (2010). 

22 D. Kim, S. Cho, N. P. Butch, P. Syers, K. Kirshcnbaum, 
S. Adam, J. Paglione, and M. S. Fuhrer, arXiv: 1105. 1410 
(2011). 



